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OBSERVATIONS OF PROTON SPECTRA (1.0 < E < 300 keV) AND 

p 

PITCH ANGLE DISTRIBUTIONS AT THE PLASMAPAUSE 
D. J. Williams, T. A. Fritz, and A. Konradi 

Detailed proton spectral and pitch angle distribution 
observations have been obtained from instrumentation flown 
aboard the NASA Small Scientific Satellite, S 3 -A (Explorer 
45). S 3 was launched from Kenya, Africa, into an elliptical 
orbit having an apogee of 5.24 R^, a perigee of 220 km, an 
inclination of 3.5°, and a period of 7.82 hours. Data used 
in this initial report are from two proton detectors and a 
three axis fluxgate magnetometer. Data from the magnetometer 
are used to routinely display the particle data as a function 
of local pitch angle. The low energy proton instrument con- 
sists of an electrostatic analyzer-channeltron configuration 
measuring proton energies from 0.8 keV to 28 keV in 16 energy 
intervals. The high energy proton detector is a telescope 
detector system consisting of two surface barrier solid state 
detectors behind a 2.2 kilogauss magnet used to sweep out 
electrons of energy less than 300 keV. Electrons of energy 
> 300 keV are detected via the coincidence mode of the tele- 
scope. This latter instrument measures protons from 24 keV 
to 300 keV in six intervals. All pre-flight and in-flight 
calibrations indicate completely nominal operation. The sat- 
ellite spin axis was placed in the orbital plane, and all 



pitch angle distributions are obtained by sectoring the spin 
cycle into a number of samples. For the data reported herein, 
32 samples per spin period (8.451 seconds) were obtained. See 
Longanecker and Hoffman (1972) for a more detailed description 
of the satellite and instrumentation. 

The data of interest are from S 3 orbit 99 in-bound occur- 
ring on December 17, 1971, some 8 hours prior to the sudden 
commencement of a magnetic storm occurring on that date. The 
magnetic conditions at the time are shown in Cahill (1972) 
where moderate magnetic activity (30-40 y) is present from 
apogee to L ^ 3 on orbit 99. Figure 1 shows proton intensi- 
ties in the energy range 24.3 to 35.1 keV as a function of 
altitude. Data are shown for two equatorial pitch angles, 

90° and 45°. The solid vertical bar at the top of the figure 
indicates the point at which the S 3 electric field detector 
comes out of saturation and gives a rough estimate as to the 
location of the high density region of the plasmapause 
(Maynard and Cauffman, 1972). Three features to note in 
figure 1 are: 

1. The high altitude (L > 4.0 earth radii) region of 
this particle distribution displays a relatively 
normal pitch angle distribution in that the inten- 
sity versus pitch angle plots peak at local pitch 
angles of 90°. The magnetic latitude of S 3 at this 
time is ^12°. The data indicate that at these 
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higher altitudes the pitch angle distributions dis- 
play a general form of sin n a with n varying from 
1/3 - 2 depending on spatial location and energy. 
Variations from this form are often seen, especially 
at small pitch angles, but will not be discussed here. 

2. At the steep radial gradient on the inner edge of 
the proton distribution the pitch angle distributions 
go over into a non-normal mode displaying strong off- 
equatorial peaking. This is evidenced by the fact 
that many more particles are observed at a pitch 
angle of 45° than at 90°. 

3. The region of the large radial gradient corresponds 
well with the plasmapause location as measured by 
the electric field detector. 

We have studied the evolution of proton pitch angle dis- 
tributions for all energies measured throughout the spatial 
regions shown in figure 1. A sample of this evolution is 
presented in figure 2. Here the proton intensity is shown as 
a function of local pitch angle for three energy bands situ- 
ated in the heart of the ring current distribution. Pitch 
angle distributions are shown from the high altitude "normal" 
portion of the particle distribution on into the gradient 
region shown in figure 1. We see that as the altitude de- 
creases, the two high energy bands lose particles at all pitch 
angles. However, particles mirroring in the near equatorial 
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Proton Intensity 24.3 Kev < E p < 35.1 Kev 
Dec 17, 1971 S 3 Orbit 99 Inbound 



20 30 40 50 

L (Earth Radii) 


Figure 1. Proton intensity as a function of altitude. Sec- 
tor = 0.264 seconds. Solid vertical bar at L = 3.45 is 
where DC electric field detector comes out of saturation t 
thereby providing indication of plasmapause (Maynard and 
Cauffman , 1972). 
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Selected Proton Pitch Angle Distributions 
S 3 Orbit 99 Inbound 



Local Pitch Angle (degrees) 

Figure 2. Evolution of proton pitch angle distributions through region of 
steep gradient s as S 3 was at geomagnetic latitude of ^12°. 




regions are lost preferentially with the result that in the 
spatial gradient region of the particle distribution, pitch 
angle distributions at the higher energies show many more 
particles mirroring at small pitch angles than in the near 
equatorial regions. The lower energy band begins to display 
the equatorial minimum in the pitch angle distribution at the 
lowest altitude shown in figure 2. Note that even though the 
center energies of the two lowest channels shown in figure 2 
are separated by only 4 keV, a sharp difference exists in the 
evolution of the pitch angle distribution as a function of 
altitude. This difference we feel is due to two causes: (1) 

a resonant energy effect and (2) strong energy diffusion. 

In figure 3, we present the differential flux values as 
a function of altitude for a = 90° ± 10°. The two highest 
differential energy bands display a radial profile which is 
qualitatively quite similar to radial profiles displayed pre- 
viously in the literature for integral energies above 'x.l 00 
keV (Davis and Williamson, 1963; Williams, 1970). However, 
note that the lower* four energies show a markedly different 
behavior. As the altitude decreases the lower four channels 
are either at or increase to a maximum flux value and at some 
lower altitude suddenly decrease in intensity resulting in a 
steep low altitude intensity gradient. It is also seen that 
the lower differential energy channels show the intensity de- 
crease at lower altitudes. This is in qualitative agreement 
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with the initiation of the ion-cyclotron instability in which 
the resonant energy is inversely proportional to the total 
plasma density and would thus decrease as the plasmapause 
(higher density) is approached (Kennel and Petschek, 1966; 
Cornwall et al . , 1970) . 

In figure 4, we present the evolution of the proton dif- 
ferential energy spectrum for a = 90° ± 10° as the region of 
the spatial gradient is approached. Three points of interest 
in figure 4 are: 

1. The development of a deep trough in the differential 
energy spectrum; 

2. the appearance of a high energy peak in the differ- 
ential energy spectrum; and 

3. a build-up of intensities observed at energies just 
below the strongly resonant region indicative of sig- 
nificant energy diffusion especially at L = 3.6 and 
from 10-20 keV. 

We wish to present the following discussion arguing that 
the preceding data are consistent with the initiation of the 
ion-cyclotron instability which causes the initial decrease in 
intensities as the region of the radial gradient is approached. 

Exponential wave growth is expected if the energetic pro- 
ton population satisfies the following conditions (Kennel and 
Petschek, 1966; Cornwall et al., 1970): 
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Proton Differential Energy Spectra 
S 3 Orbit 99 Inbound 



Figure 4. Evolution of proton differential energy speatra 
through region of steep gradient } a = 90° ±10°. 
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A = a measure of proton pitch angle anisotropy 
(A s ^ for sin n a distribution) 
ft = proton gyrofrequency 
to = wave frequency 


(2) J (E > E D ) > KL" 4 protons/cm 2 sec 

Max K 

L = Mcllwain L parameter 
E r = parallel resonant energy 


(3) 
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or from ( 1 ) 

e r 5 §¥n a ' 2 (1+A)-‘ 

B = magnetic field magnitude 

N = total plasma density 

Similar criteria hold for the stimulation of such insta- 
bilities in the energetic electron population (Kennel and 
Petschek, 1966) and the effects of the electron branch of the 
instability (whistler mode interaction) have been discussed 
in connection with electron precipitation phenomena (Brice, 
1970; Brice and Lucas, 1971). 
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The pitch angle anisotropy identifies at a given point 
in space the wave frequencies which will resonate with the 
proton particle distribution. We thus interpret the initia- 
tion of the intensity decreases in figure 3 in the following 
manner. The lower four differential channels show proton 
fluxes at high altitudes which exceed the limiting flux given 
in (2). As the altitude decreases, the resonant energy de- 
creases due to an increasing plasma density as the plasmapause 
is approached. When the resonant energy matches a particular 

differential energy band, protons in that band suffer an in- 
tensity decrease due to the initiation of the ion-cyclotron 
instability as all conditions are then met. As the altitude 
continues to decrease (N increasing) this effect is progres- 
sively seen through the lower energy channels as displayed 
in figure 3. The high energy channels do not show similar 

effects because they do not exceed J u 

Max 

To more quantitatively test this hypothesis, we show in 
figure 5 omnidirectional flux values versus altitude. Here 
we have summed through the lowest four energy channels which 
appear to be in resonance and have shown the higher energy 
channel by itself which is apparently unaffected. For com- 
parison, we have included a plot of the maximum allowed flux 
value from (2) for two values of the constant K, 5 x 10 10 and 
10 10 cm" 2 sec _1 . From figure 5 we see that a critical flux 
value of 5 x 10 10 L -4 protons/cm 2 sec , is consistent with the 
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J (Protons/cm 2 sec) 



L (Earth Radii) 


Figure 5. Omnidirectional fluxes versus altitude compared with J max = KL ~ 4 fo 
two values of K . Vertical bar at L = 3.45 3 same as in figure 4, a = 90° ±10 



above discussion. Theoretical estimates of K yield a lower 
limit of 5-7(10) 10 cm _2 sec _1 (Kennel and Petschek, 1966; 
Cornwall, 1966) and recent low altitude experimental normali- 
zation of K has given 'v 1 . 3 ( 1 0) 1 °cm~ 2 sec~ 1 (Cornwall et al., 
1971a). 

Utilizing this interpretation we then are able to measure 
the resonant energy within a certain differential energy band 
as a function of altitude from figure 3 by obtaining the 
altitude at which the sudden intensity decrease occurs. We 
therefore can extract the total plasma density N, as a func- 
tion of altitude by observations of the resonant energy, E D , 
the pitch angle anisotropy. A, and the magnetic field, B. 

This result is shown in figure 6. The vertical bars are the 
values of N extracted and reflect the fact that the resonant 
energy is measured as a differential energy band. A dashed 
line is drawn through the vertical bars to indicate that a 
normal plasmapause structure is obtained from the interpreta- 
tion that the data in figures 3, 4, and 5 are consistent with 
the i ni ti ati on of the ion cyclotron instability. 

It is possible to explain qualitatively the evolution of 
pitch angle structures shown in figure 2 by the use of a 
resonant instability provided that energy diffusion is signif- 
icant. Excluding the singular point at a Q = , there are 

several reasons why the interaction process will be more ef- 
fective in the near equatorial regions (Kennel and Petschek, 
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Figure 6 . Plasma densities obtained from resonant energy equation with Er B a 
A measured at or gust before point of intensity decrease . N is extracted from 
equation shown in figure . Vertical range of N due to observation of Er in par 
txcular energy interval. Dotted line simply indicates that a "normal" plasma- 
pause xs consistent with the data. 




1966). For example, wave growth rate is largest near the 
equator due to the small B and thus lower resonant energies 
required. In addition, particles mirroring well off the 
equator spend less time in the interaction region (have large 
v 1 1 at the equator) and are thus less effectively scattered. 
Thus while particles are lost from the detecting channel 
through both pitch angle and energy diffusion at all pitch 
angles the effect is greater near the equator. 

If the above interpretation is correct and protons are 
put into strong pitch angle diffusion, the intensities shown 
in figures 3 and 4 represent sufficient energy to cause a 
mild red arc through the mechanism suggested by Cornwall et 
al. (1971b). In checking the records of the Fritz Peak Ob- 
servatory of the NOAA Aeronomy Laboratory, it was found that 
a mild red arc was observed at the time these observations 
were made (Hernandez, personal communication). Observations 
of the red arc occurred at near local midnight and the radial 
gradient displayed in figure 1 was observed as S 3 passed 
through the local midnight meridian. Rough estimates have 
shown that the energy available to the ionosphere from the 
particles shown in the preceding figures was sufficient to 
support the intensity of the red arc observed. 

For the data presented in this report, w/fl. is approxi- 
mately 1/5. At L = 4, this gives an Alfven wave frequency 
cutoff of approximately 1.6 Hz. In the gradient regions. 
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3.4 < L < 3.6, enhanced variations in the 1 to 3 Hz channel 
of the electric field detector were observed (Maynard and 
Cauffman, 1972). However, no variations were seen above the 
sensitivity threshold for the 1 to 3 Hz channel of the search 
coil magnetometer (Parady, personal communication). It is 
possible that with the intensities involved, the growth rate 
would be such that wave amplitudes exceeding the search coil 
magnetometer thresholds would not have occurred. 

In summary, we wish to point out that the data presented 
here are consistent with the initiation of the ion cyclotron 
instability when the requirements (1), (2), and (3) are met. 
These requirements are met initially at the altitude at which 
the sudden intensity decrease occurs (figure 3). However, 
after the initiation of the instability, we find that the 
linear theory is outlined by (1), (2), and (3) is unable to 
explain the further evolution of the intensities, pitch angle 
distributions and energy spectra of the ring current particles. 
First, the theory as it now stands has no mechanism whereby 
energy loss is formally introduced into the scattering mech- 
anism. Thus the pitch angle distribution of figure 2 cannot 
be satisfactorily explained via ion-cyclotron resonance. 
Secondly, the present theory also has no mechanism whereby 
the fluxes would precipitate much below J . The instability 
should self-quench and basically return the particle fluxes 
to or just below J max * However, we see from figures 3, 4, 
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and 5 that this is not the case and the intensities continue 
to decrease to one to two orders of magnitude below this 
critical flux limit. Therefore to describe the actual evolu- 
tion of the data presented herein, we feel that non-linear 
effects, additional instabilities, and particle injection and 
drift as a function of pitch angle (Konradi et al . , 1972) 
must be considered in more detail. 
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